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Abstract 

We explore the statistical behavior of the order statistics of the flights 
of One-sided Levy Processes (OLPs). We begin with the study of the 
extreme flights of general OLPs, and then focus on the class of selfsimilar 
processes, investigating the following issues: (i) the inner hierarchy of the 
extreme flights - for example: how big is the 7*'' largest flight relative 
to the 2""^ largest one?; and, (ii) the relative contribution of the extreme 
flights to the entire 'flight aggregate' - for example: how big is the 
largest flight relative to the OLP's value?. Furthermore, we show that 
all 'hierarchical' results obtained - but not the 'aggregate' results - are 
explicitly extendable to the class of OLPs with arbitrary power-law flight 
tails (which is far larger than the selfsimilar class). 

Keywords: Levy flights, selfsimilar Levy processes, order statistics, 
extreme value theory, Frechet distribution. 

*This manuscript is an extended version of a contribution to a special Physica A volume 
in honor of Shlomo Havlin on his sixtieth birthday. 
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1 Introduction 



Levy processes - random motions with stationary and independent incre- 
ments - constitute one of the most important and fundamental family of stochas- 
tic processes. Special examples of the Levy family include Brownian motion 
(Wiener process), Poisson processes, and Compound Poisson processes. Since 
their introduction in the 1930s, by the French mathematician Paul Levy pQ- 
P], Levy processes were studied extensively by both theoreticians and applied 
scientists. The literature on Levy processes is vast, and their range of appli- 
cations encompasses numerous fields of science and engineering. See |1]-|12| 
for the theory of Levy processes, and JJ-QT] and references therein for their 
apphcations. 

Amongst the family of Levy processes, the class of One-sided Levy Processes 
(OLPs), also referred to as Levy Subordinators, is of special importance. OLPs 
satisfy the additional requirement of non-negativity of their increments, render- 
ing them monotone non-decreasing and non-negative valued. OLPs are natural 
models for random flow of positive valued quantities. Examples include: mass, 
energy, and time in physical systems; work and costumers in queueing systems; 
data in communication systems; claims in insurance; etc. 

Unlike Brownian motion - whose sample paths are continuous, the sample 
paths of OLPs - regardless of their statistics - are always purely discontinuous. 
The propagation of OLPs is conducted only by flights (jumps), and not by any 
sort of continuous motion. This is due to the fact that OLPs are continuum 
superpositions, or aggregates, of Poisson processes. 

Since the structure of OLPs is that of a 'flight aggregate', questions of the 
following type arise naturally: how big is the largest flight? how big is the n*"^ 
largest flight? how big is the m"^ largest flight relative to the largest flight 
(m > n)? how big is the largest flight relative to the entire aggregate? how 
big is the combined contribution of the rt*'' largest flights relative to the entire 
aggregate? - that is, questions regarding the order statistics of the flights of 
OLPs. 

The study of the order statistics of series of Independent and Identically 
Distributed (IID) random variables is a well established field of probability 
called Extreme Value Theory (EVT). This theory originated, in the 1920s and 
30s, with the pioneering works of von Bortkiewicz ,18) . Frechet J_9,, Fisher 
& Tippett 1201, von Mises Ell, WeibuU, and Gunibel ^ (the theory's three 
possible types of limiting distributions are named after Gumbel ("type 1"), 
Frechet ("type 2"), and WeibuU ("type 3")). A rigorous theoretical framework 
was presented in 1943 by Gnedenko [23]. The statistical analysis of extreme 
values is of major importance in the analysis of rare and 'catastrophic' events 
such as floods in hydrology, large claims in insurance, crashes in finance, material 
failure in corrosion analysis, etc. See [23" EH] Etnd references therein for both the 
theory and applications of modern EVT. See also [2Z] for a recent application 
to the study of complex networks. 

OLPs are the continuous time counterparts of non-negative IID random 
sequences (viewed as discrete time stochastic processes). When passing from 
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discrete time to continuous time the question "what was the largest observation 
amongst the first N observations?" transforms to "what was the largest flight 
that occurred up to time T?" . More generally, the analogue of the order statis- 
tics of a discrete sequence (largest observation, second largest observation, etc) 
is the order statistics of the flights of a continuous time process (largest flight, 
second largest flight, etc). Hence, the study of the extreme flights of OLPs is 
also a natural sequel to the EVT of IID random variables. 

Furthermore, the only possible non-trivial scaling limits of non-negative IID 
random sequences are the, so called, selfsimilar OLPs. This special class of mo- 
tions, which occupies a predominant role in both the theory and applications of 
stochastic processes, is composed of all OLPs which are invariant under changes 
of scale. Since EVT is, in essence, an asymptotic theory (providing limiting 
results as the total number of observations N tends to infinity) we shall devote 
our main emphasis in this work to the exploration of this special class of OLPs. 

The paper is organized as follows; We begin, in section [3 with an investi- 
gation of the extreme flights of general OLPs (following a short review of these 
processes). In sections|31and0]we focus on the class of selfsimilar OLPs. Section 
|3|studies of the internal hierarchy of the extreme flights - dealing with questions 
of the type "how big are the extreme flights relative to each other?" . Section 

01 on the other hand, studies the magnitude of the extreme flights relative to 
the entire process - the ^flight aggregate\ Furthermore, we explain why (and 
how) the results of section O are extendable to the entire class of OLPs with 
power-law flight tails, and, on the other hand, why the results of section 01 do 
not extend. 

Throughout the manuscript: P(-) = Probability; and, E[-] = Expectation. 

2 OLPs and their extreme flights 

In this section we give a short review of OLPs, establish a few preliminary results 
regarding the distribution of their extreme flights, and conclude with the special 
class of selfsimilar OLPs and its connection to the Frechet ("type 2") extreme 
value distribution 25 . 

2.1 Levy processes and OLPs 

A stochastic process is said to be Levy if it is continuous in probability and 
has stationary and independent increments. A OLP is a Levy process with 
non-negative increments. 

Characterization of OLPs 

The celebrated Levy-Khinchin formula (see, for example, jl] or |5]) asserts 
that Levy processes can be decomposed into two independent parts: (i) a purely 
continuous part, which is a Brownian motion; and, (ii) a purely discontinuous 
part, which is a superposition of Poisson processes. Since Brownian motion 
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is symmetric, OLPs can not include a continuous part and are hence purely 
discontinuous (pure jump) processes. 

A OLP L = (i(i))(>o is characterized by its Laplace transform (its 'spectral 
representation' in Laplace space): 

E[exp{-wL(t)}] =exp{-^'(w)-t} ; w>0. (1) 

The function ^' is called the Levy characteristic of the OLP L (in the literature, 
^ is also referred to as the spectral characteristic, or symbol, of L). 

The Levy measure 

If L is a Poisson process with flights of size xq {xq > 0) and rate Aq, then 
its Levy characteristic is ^(w) = (1 — exp{— wxo}) Aq. If L is a superposition 
of N independent Poisson processes - process n, n — 1,2, ■ ■ ■ , N , having flights 
of [xn > 0) and rate A„ - then its Levy characteristic is given by 

N 

4'(w) = ^ (1 -exp{-wx„})A„ . (2) 

n=l 

Hence, passing from (O to a continuum limit, where flights of size x {x > 0) 
occur at rate X{dx), we arrive at 

/"OO 

*(w) / (1 -exp{-a;a;})A(dx) . (3) 
Jo 

Equation (^J is the Levy-Khinchin representations (in Laplace space) for 
OLPs. The rate A(-) - the Levy measure of i - is a measure on the non- negative 
half line (0, oo) satisfying the integrability condition min{a;, l}A(da;). The 
total rate X{dx) could certainly be infinite - not due to non-integrability at 
X — oo but, rather, due to possible non-integrability at a; = (intuitively, large 
flights can occur only rarely, but tiny flights may occur very frequently) . A OLP 
is Compound Poisson if and only if its Levy measure has finite total mass. 

Given a Levy measure A(-) it is natural to introduce its 'cumulative distri- 
bution function'. However, since A(-) might have inflnite total mass (due to 
possible non-integrability at x = 0) we need to define the 'cumulative distribu- 
tion function' by integrating (backwards) from x — oo rather than (forward) 
from a; = 0: 

/•oo 

A{x) = / X{dy) ; X > . (4) 

J X 

We henceforth refer to A as the OLP's flight tail. The meaning of the flight tail 
is straightforward: A{x) is the rate at which flights of size > x occur. 

Examples 

We give a few examples of classes of OLPs. The first two have finite Levy 
measure, whereas the other have infinite Levy measure; 
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1. Compound Poisson with exponentially-distributed jumps (a > 0): 



X{dx) = aexp{—ax}dx (A(a;) = exp{— aa;}) , 
^ ^ a + (j 

2. Compound Poisson with Gamma-distributed jumps {a,p > 0): 

X{dx) = exp{—ax}xP~^dx , 

^^!{io) = ^{p)(—-^^—-] . 

3. Gamma processes - OLPs with Gamma-distributed increments (a > 0): 

\iA \ exp{-aa;} 
Alaa;) = -dx , 

X 

= In (^1 + ^) . 

4. Selfsimilar ('fractal') OLPs (0 < a < 1): 



a / 1 



X" 



^{oj) = r(l - a)w°' . 
2.2 The extreme flights of OLPs 

Given a OLP L = {L{t))^yg with Levy measure A(-) (flight tail A(-)), we set 

Xi{t) > X2{t) > Xait) > ■ ■ ■ (5) 

to be the order sequence of the flights of the OLP L, i.e; X„(<) denotes the n*^^ 
largest flight that occurred during the time interval [0,t]. If the OLP has finite 
Levy measure A(0) < oo (i.e; if it is Compound Poisson) then, up to time t, 
there will have occurred only a finite number of flights (Poisson with rate tA(0), 
to be exact) and the order sequence will hence be finite. On the other hand, if 
the OLP has infinite Levy measure A(0) = oo then the order sequence will be 
infinite. 

Since the Poisson distribution will appear time and again in the sequel, we 
introduce the following shorthand notation for its probability frequencies: 

Pn{t^) = ^eM-t^) ■ (6) 

That is; Pn{l^) denotes the probability that a Poisson random variable with rate 
l-i {n > 0) equals the integer value n (n = 0, 1, 2, • • • ). 
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The joint distribution of the flights of a OLP is of major importance. Given 
a time t {t > 0) and an interval / C (0,oo), let n(t; /) denote the number 
of flights of L, during the time period [0,t], whose size laid in the interval /. 
Due to the 'Poissonian-superposition' structure of L the random variable Il{t; I) 
is Poisson-distributed with rate t Jj X {dx) . Moreover, if /i , ■ • ■ ,Ik are disjoint 
intervals then n(i; /i), ■ • ■ , n(i; Ik) are independent random variables and hence 



P (n(i; Ik)=nk;k^l,--- ,K)^Y[ P„, (t[x {d: 

fc=i ^ 



(7) 



Vni, • • • , UK E {0, 1, 2, • • • }. For further details see 
The greatest flight 

The distribution of the maximal flight Xi{t) is given by: 

P {Xi{t) < = exp{-tA{x)} . (8) 
The explanation of (|SJ) is given by the following straightforward deduction: 

P{Xi{t)<x) 

= P(n(i;(x,^)) = 0) 

= Po{tCXids)) 

= exp{-<A(a::)} . 

Equation © enables the 'shock tolerance' design of systems. We explain; 

Assume that L = (L(t))j>Q is the inflow to a system, which cannot absorb 
flights (i.e; inflow surges) of size greater than a tolerance level I. The following 
natural engineering question arises: what should the tolerance level I be so that 
to ensure that the system would withstand the time period [0, T] with no failure, 
with probability greater than 1 — S [0 < S < 1)? The answer, deduced by a 
simple inversion of ||SJ) is 

/^A-rih/ ' 



T 

The maximal flight, viewed as a stochastic process {Xi{t))^yQ, is Markovian. 
Indeed, when at level x {x > 0): (i) the process has to wait an exponential 
time, with rate A(a;), till transition; and then, (ii) the process will transit to the 
level y {y > x) with probability A (dy) /A{x). In other words, the inflnitesimal 
generator of the process (Xi(t))j>Q is given by the integral operator: 

(Afyf) (a;) = / {(p{y) ~ (fix)) X{dy) . 



The n runner-up 
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The distribution of the n^^ largest-order flight X„ (t) is given by: 

P {Xn{t) < a;) - ^ Pfc {tA{x)) . (9) 

k=Q 

The explanation of is analogous to the deduction of 0: 

P{Xr,{t)<x) 

= P(n(t; (a;,oo)) < n - 1) 

The order sequence 

The entire order sequence can be constructed/simulated sequentially, ac- 
cording to the iterative scheme (y < x): 

P {Xn+i{t) < y I X„(t) = x) - cxp{-t (A(y) - A{x))} , (10) 

together with the initial condition ||SJ|. Indeed, Vy < x we have 

P {Xn+l{t) < y \ Xn{t) = X) 

= P(n(t; (y,x)) =0) 

= exp{-t(A(y)-A(a;))} . 

Furthermore, the joint distribution of the truncated order sequence {Xi{t), • • ■ , 
is given, for a decreasing sequence oo > 6i > ai > 62 > a2 > ■ • • > > qk = 0, 
by the following formula: 

P{ak <Xk{t) <bk ; fc = I,-- - ,X) 

(11) 

eM-tA{bK)} ■ t^-' UkJi (A(afc) - A{b,)) 
The derivation of (|ll|l is analogous to the derivation of l|5|l- (|10|) . 

2.3 Self-similarity and the Prechet distribution 

In the proceeding sections we shall focus on selfsimilar OLPs, i.e; OLPs with 
flight tail A{x) ~ l/x" (0 < a < 1). For these processes the distribution of the 
maximal flight Xi{t) is given by (due to ||H1)): 

P {Xi{t) <x)= exp{-t/x"} . (12) 
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The Prechet distribution 

The probability law H12|l is known as the Frechet distribution, or as the "type 
2" extreme value distribution It emerges as the asymptotic probability law 
of the maximum of Independent and Identically Distributed (IID) heavy tailed 
random variables. We explain; 

Let Yi, - ■ ■ , Y/v be a sequence of IID, non-negative, random variables with 
heavy ('fat') tails: P{Y > y) ^ a/y"' as j/ ^ cx) (for some a,a > 0). Then, 
the appropriately scaled sequence maximum converges, in law, to the Frechet 
distribution: 

^ max{yi, ■■■,YN}<y^ cxp{-a/y"} . (13) 

Indeed; 

P(^max{yi,--- ,YN}<y) 

= (1 - p (r > Tvi/"?/))^ 



N 



~ exp{-a/y"} . 

Note, however, the discrepancy between (|12|l and H13|) : in H12|) we are re- 
stricted to the parameter range^ < a < 1, whereas in (|13|l the exponent a may 
admit any positive value. This discrepancy stems from the following reason; 

Consider the IID sequence Yi, ■ • ■ ,Yn introduced above. The scaling limit 
of its maximum is given by 1)13(1 . for all a > 0. However, in the scaling limit of 
their sum (aggregate) Sn = Yi + - ■ ■+Ym a phase transition occurs when passing 
from a < 1 to a > 1: (i) for a < 1 the appropriate scaling is Sn/N^^" and 
the limiting distribution is a selfsimilar Levy law of order a; but, (ii) for a > 1 
the scaling is the universal Law of Large Numbers scaling - Sn /N - leading to 
the deterministic limit E[y]. Hence, scaling limits of non- negative IID sums - 
yielding the selfsimilar OLFs - are restricted to the 'Levy range' < a < 1, 
whereas scaling limits of non-negative IID maxima admit the entire 'Frechet 
range' a > 0. 

This gap between the 'Levy range' and the 'Frechet range' can be bridged 
by taking scaling limits of OLPs with power-law flight tails; 

Scaling limits 

Let i be a OLPs whose flight tail asymptotics (as x — > oo) are A(x) ^ l/a;"^, 
for some exponent a > 0. Introduce the scaled order sequence, x''i\t) > 
X^'\t)>X^^'\t)>---, given by 



^Selfsimilar OSLMs with exponent a > 1 are trivial: L{t) = const -t if a = 1, and L{t) = 
if a > 1. 
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(where c > is a scaling factor), and consider the scaling limit of the order 
sequence, i.e; the hmit, in law, as c ^ oo, of the scaled sequence: x'f"\t) > 
X^^\t) > X^°°\t) > ■ ■ For the scahng hmit x[°°\t) we have: 



(^X[°^'\t) < x'^ =exp{-i/a;"} 



and now the exponent a is not (!) restricted to the 'Levy range' < a < 1 - as 
it is in ((T^ . 

In fact, we obtain the following counterpart of Hll|) : 



(14) 



0° (bO 



which holds for the entire 'Frechet range' a > 0. The deduction of p4|) stems 
from combined with the scaling procedure described above. Indeed, for 
c >> 1 we have: 



P (^ak < Xj^\t) < bk ; fc = l,-- - 

= P {cak < X„(c"i) <cbk ; k^l,--- ,K) 

= eM-{c"t)A{cbK)} ■ (c"t)^-i Ylk=i (A(cafe) - A(c6fe)) 



3 The hierarchy of extreme flights 

In this section we explore the hierarchical structure of the order sequence of 

selfsimilar OLPs. We assume, throughout this section, that the OLP L is a- 

selfsimilar, i.e; that its flight tail is A(x) = l/x". 

The Beta distribution will emerge naturally and play a key role in this section 
We denote by Bn^m{x), x G [0, 1], the cumulative distribution fmiction of a 

Beta{n,m) distribution {n,m > 0): 

BnM^) = r Il!i±/^«"-i(l - ur-'du . (15) 

Jo r(n)r(m) 

Let us begin with the investigation of the pairwise hierarchy; 
3.1 Pairwise hierarchy 

In this subsection we study the statistics of the ratio Xn+Tn{t)/Xn{t), i.e; the 
size of the {n + m)*'' largest-order flight relative to the size of the n^^ largest- 
order flight. Clearly, this ratio admits values ranging in the unit interval. We 
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assert that the cumulative distribution function of this ratio is independent of 
time t and is given,Vu G [0, 1], by: 



(16) 



Equation H16|) is a particular case of proposition ^ which will be presented in 
the following subsection. 

Moreover, the moments of the ratio Xn+m{t) / Xn{t) are derived from (|16|) : 
Vp > we have: 



E 



X, 



7i-\-7n 



n+m — 1 

n 



k—n 



k + p/a 



The proof of IjlTfl is brought in the appendix. 
We point out two special cases: 

'Consecutive relativity' 

The statistics of the ratio Xn+i{t) / Xn{t) are given by: 



and 



Xn{t) 



< U 



-p/c 



(17) 



'MsLximal relativity' 

The statistics of the ratio Xi+„i(i)/Xi(t) are given by: 

Xl+m{t) 



Xi{t) 



<u] =1- (1 -u") 



and 



E 



X 



1+m 



Xi{t) 



n 

fe=i 



k 



k + p/a 



3.2 Multidimensional hierarchy 

In the previous subsection we investigated the pairwise hierarchy of the extreme 
flights. Now we turn to study the multidimensional hierarchical structure of the 
order sequence: given an increasing sequence of integers ni < ■ • • < uk < nx+i 
we wish to compute the joint distribution of the vector of ratios 



Xn^iit) Xn^{t) Xn^^i{t) 



Xni (t) Xn2 (t) 



(18) 



To that end we have; 
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Proposition 1 The multidimensional cumulative distribution function of the 
vector of ratios ^18\) is independent of time t and is given, Vui, • • • , uk G [0, 1], 
by 

The proof of proposition^ is brought in the appendix. This proposition can 
be re-statcd as follows: 

The ratios in jJij^) are independent and are distributed according to fldl) . 



In particular, proposition ^ implies that: 
Xk{t) 



P ( < «fe ; fc = 1, • • • , if) = n iu,r' , (20) 

^ k=l 



VUI, • • ■ G [0, 1]. 

Multidimensional 'maximal relativity' 

The distribution of the truncated order sequence relative to the maximal 
flight 

'X2{t) x^it) XK+iity 



,xi(t)'xi(t)' ' xi(t) ; ' ^^^^ 

rather than its multidimensional 'consecutive relativity' counterpart (|2U|I . is 
given, for a decreasing sequence 1 > oi > • • • > ok > ok+i = 0, by the 
formula: 



-xT{t) 

The proof of H22|) is brought in the appendix. 



(22) 



A 'Prechet remark' 

All the results of this section may be extended from the restricted selfsim- 
ilar 'Levy range' < a < 1 to the entire 'Frechet range' a > following 
the scaling procedure described in subsection 12.31 replace the order sequence 
Xi{t) > X2{t) > Xsit) > ■■■ oi an a-selfsimilar OLP (0 < a < 1) by the 
scahng limit x[°°\t) > X^°^\t) > X^^°°\t) > ■■ of the order sequence of a 
OLP whose flight tail asymptotics are A(a;) ~ l/x" (a > 0). 
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4 The extreme flights vs the aggregate 



In the previous section we studied the behavior of the extreme flights relative 
to each other. We now turn to study the extreme flights relative to the entire 
process L. Since L is a (continuum) superposition of Poisson processes, it is, in 
fact, a ^flight aggregate^ (where flights of size x occur with rate X(dx)). Hence, 
we are interested in investigating the relative contribution of the largest flights, 
up to time t, to the entire 'fight aggregate' L{t). As in the previous section, we 
assume in this section that the OLP L is a-selfsimilar. 



Before we begin, let us introduce the family of functions {Gq}q^^^-^, defined 
on the non- negative half line (w > 0), and given by the integral formula 

Ga(w) = a / ^ ^-du , (23) 

Jo " 

or, equivalently, by the power series 

G.(..) = -a^Li^.i^. (24) 

m—l 



4.1 The greatest flight vs the aggregate 

Let ^{t) denote the magnitude of the "entire flight aggregate, up to time t, 
except for the maximal flight", relative to the maximal flight: 

_ L{t)-Xi{t) 

x,{t) ■ 

We assert that the ratio ^{t) is independent of time t and that its Laplace 
transform is given by: 

E[exp{-^g(t)}]= 1, . (25) 

J- + '-Jay^^j 

Equation (|25|l is a particular case of proposition|21which will be presented in the 
following subsection. Let us denote by Qa the probability law whose Laplace 
transform is given by the right hand side of (|25fl . 

Equation H25|) implies that the mean and variance of ^ {t) are independent of 
time t, and are given, respectively, by 

E[^(i)] = -^, (26) 
1 — a 



and 



^-(^(^»= (l-.)^(2-a) - 
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The derivation of H26|) - (|27|) is obtained by differentiating H25II and using the 
power expansion 1)241) . Combining (|26() - (|27|l together we also have 

Var ( J}^] = ,J , . (28) 

From ||in|)-(|lHll we see that: 

When a ^ 1 the underlying OLP L converges to the degenerate determin- 
istic linear limit = t, and hence all flights tend to zero. Indeed, from H26|l we 
obtain that E[^(t)] ^ oo as a ^ 1. That is, the size of the aggregate, relative 
to the greatest flight, tends to infinity. Moreover, the variance of normalized 
ratio ^(t)/E[^(i)] converges, as a ^ 1, to 1. 

On the other hand, when a ^ the tails of the OLP L become infinitely 
'fat' and hence the greatest flight should dominate. And indeed, we obtain that 
E'K(i)] — > as a — > 0. That is, the size of the aggregate, relative to the greatest 
flight, tends to zero. In fact, the dominance of the largest flight is so great that 
even Var(^(t)) — > as a 0. 

The value a = 1/2 turns out to be the "break-even" point where the mean 
of ratio ^(i) equals 1. That is, the entire aggregate splits, on average, half-half 
between the largest flight and "all the rest". As a 'moves up' towards 1 the 
weight of "all the rest" takes the lead and eventually dominates, and vice- versa 
as a 'moves down' towards 0. 

Another immediate consequences of (|25|) . after noting that Xi{t)/L{t) = 
1/(1 -I- ^(t)), is that the cumulative distribution function of the ratio Xi{t) / Lit) 
is independent of time t and is given, Vu G [0, 1], by 

■'^■'"<.)^pf«>i-ll . (29) 



L{t) ~ ) \ u 



where ^ is ^^-distributed. 



4.2 The top-n vs the aggregate 

We now proceed to study the contribution of the 'top-n' flights, Xi{t),X2{t), ■ ■ ■ ,X„ 
relative to the entire aggregate L{t). The key result is: 



Proposition 2 The random vector 

' L{t) ~ {Xi{t) + ■ ■ ■ + X„+i{t)) Xi{t) XAt) 



(30) 



Xn+l{i) ' Xn+l{t)^ ^ Xn+l{t) 

is independent of time t and equal, in law, to the vector 

where: (i) S^^, • • • , Yi, • • • ,Yn are independent random variables; (ii) Cii ' ' ' i 

are Q a- distributed; and, (Hi) ^(i), • • • , ^(n) are the order statistics o/ Yi, • • • , y„ 
which, in turn, are Pareto{a) -distributed: P{Y > y) = l/y", 2/ > 1- 
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The proof of proposition |21 is brought in the appendix. 



Noting that Y/^i-j + • ■ • + Y(„) ~ Yi + • • • + y„, proposition [21 implies the 
following pair of corollaries: 

(i) The ratio Xn+i{t) / L{t) - the contribution of the {n + 1)**^ largest flight, 
relative to the entire aggregate - is equal, in law, to 

{i + ^i + --- + ^n+i + Yi + --- + Y^y' . 

Hence, Vu £ [0, 1] we have: 

^^f^<u^ ^p(^^, + --- + ^,,+,+Y, + --- + Y^>^-?j . (31) 

(ii) The ratio {Xi{t) + • • • + Xn(t)) /L(t) - the combined contribution of the 
n largest flights, relative to the entire aggregate - is equal, in law, to: 

^ i + ei + ---+e„+i y' 

Yi + --- + Y„ J 
Hence, Vu G [0, 1] we have: 

,^X,(t) + ...+X„(t)^ X^p/l + er + --+e„,.-,l_^A (32) 



m -J V Yi + ---+Y,, 



A Prechet remark 

The results of this section are not extendable from the restricted selfsimilar 
'Levy range' < a < 1 to the entire 'Frechet range' a > 0. This is since the 
scaling procedure described in subsection l2.3l - which holds for any a > when 
regarding the extremes - fails to hold for the aggregate when a > 1 (since, as 
explained in subsection 12.31 the limiting behavior of the aggregate undergoes a 
phase transition when passing from the 'Levy kingdom' a < 1 to the 'Law of 
Large Numbers kingdom' a > 1). 

5 Appendix: proofs 

5.1 Equation (jl7|) : moments of the pairwise hierarchy 

Proof. Combining H15|l and (|16|l . the probability density function of the ratio 
Xn+mit)/Xnit) is givcu by 

1 (n 1 m 
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Hence, changing the integration variable to x = and using Beta integrals, we 
obtain 



E 



Jo uP ■ fn,m{u)du 

fl „(n+p/Q)-l/-| _ ^^m-l J„ 

r(n)r(m) Jo -^J "-^ 

r(n+m) T{n+p/a) 
r(n) T(n+m+p/a) 



n 



n+m — 1 



fc=n fe+p/a 



5.2 proposition [T] 

We begin with a lemma; 

Lemma 3 Let K{x) = l/x", x > 0, and n, m > 0. Then, Vx > anii VO < 
u < 1 we have 



1 



r(n)r(m) 



A(y)"-^ (A(x) - A(y))'"-i A(dy) = 



A{x 



in+m — 1 



r(n + m) 



(33) 



Proof. First, note that we can re- write the left hand side of H33() as follows 

A(x)"+™-i r fHy)V (. _ A(y) V""' hM. 

\ m) A(x) ■ 



(34) 



Now, using the change of variable s = A{y)/A{x), the integral part of H34|) equals 

/" (1 - sr-' ds = r(")r(m) ^ ^ 

Jo r(n + m) 

Substituting (|35|l back into (|34|l concludes the proof ■ 



We are ready now to prove the proposition; 
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Proof. First, we use JT)) in order to compute the multidimensional density 
function of the vector (X„j (t), - ■ ■ , (t), Xn^+i (t)) ■ 

PiXn,{t)^dxk ■,k^l,--- ,K+1) 

= ntY Pi mdxk)) Pn.^„._,-i {tQ-^ A ids)) 

= nf=Y (*A(dx.)) ( '^^^^--^-p^tL^;-"!:;""^-^"^ exp {-t (A(:.,) - A(x._i))}) 

= / (xi, • • • , a;if, a;K+i) • X{dxi) \{dxK) ■ \{dxK+i) ■ 

Vdecreasing sequence oo = xq > xi > ■ ■ ■ > xk > xk+i > 0, and where no = 0. 

Hence 

rOO \ / I \ poo \ / 7 \ 

(36) 

Jxk+1 ^ ^+1' JxK=x:K+l/y.K ^ ^' 
■ ■ ■ SZ=x,/u., ^idx2) \Z=x^lu, HdXl)f (Xi, • • • , XK,XK+l) ■ 

Now, computing the multiple integral (|36|l . using lemma |3| repeatedly, yields 



^ fe=i 



where 



1= / ^T7 rC^v{-tA{xK+l)]HxK+lT''+'~^\(dxK+l) , 

Jo T(nK+ij 
which, in turn, must equal 1. ■ 

5.3 Equation (1221) : multidimensional 'maximal relativity' 

Proof. (Recah that A{x) = 1/a;") 

Analogously to the proof of proposition ^ we have 



K 



X{dxi) J^^^^ \(dx2) ■ ■ ■ £^^1''^ \(dxK) Jq"''' KdxK+i)f (a^i, ■ ■ ■ , xk, xk+i) , 

(37) 
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where 

/ (xi, • • • , xk,xk+i) = t^~^^ exp {-tA{xK+i)} ■ 
Now, Vfc = 2, ■ • ■ , K we have 

= A{akXi) - A{ak^ixi) (38) 

= (A(afe) - A(afe_i))A(a;i) . 
And, using the change of variable s = A{xk+i), 

/p"^"^' X{dxK+i)te:Kp{-tA{xK+i)} 

= /^a^)A(xO i exp{-is}ds (39) 
= cxp{-tA(aK)A(a;i)} . 

Substituting H38|l and 13911 into the muhiple integral on the right hand side 
of ^ yields 

K 

I ■ n (A(afc) - A(afc_i)) , (40) 

fc=2 

where 

roo 

I = t^ eyip{-tA{aK)A{xi)]A{xi)^-^\{dxi) . 







However, using the change of variable s = A(a;i) gives 

I^t"" r eM-tHaK)s}s''-'ds = -^j^ (41) 
Jo A{aK)'^ 

Finally, combining H4(J|I and H41|l together implies that the right hand side 
of (jSZl equals 

^(^^ n(A(afc)-A(a,_i))= r(X)(aK)"^ n ^ ^ 



^(°^) fc=2 fc=lV(«fc+l)" (ofc)" 

■ 

5.4 proposition El 

We use the shorthand notation 

' L{t) - {Xi{t) + ■ ■ ■ + Xn+iit)) . Xi{t) Xn{t) 



iuity,vit)) 



Ar„+i(t) ' Ar„_i-i(t)' ' Ar„4.i(i) 

Proof. (Recah that X{dx) = ax^^'^+'^Ux and A(a;) = l/x") 
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step 1 

Due to the 'Poissonian-superposition' structure of the OLP L, given the 
event {Xn+i{t) = x} the vector {U{t); V{t)) is equal, in law, to the vector: 

... 

XX X 

where: 

(i) S^{t) is the value, at time t, of a OLP with Levy measure 

\{dy) < y < X 

y > X 



V{dy) = 



(42) 



(ii) J^^y • • • , J^^^ are the order statistics of Jf , • • • , J,^ which, in turn, are 
IID random variables distributed according to the probability law (supported 
on the half line {x,oo)): 



P [r G dy) 



<y < X 

X{dy)/K{x) y>x 



(43) 



(iii) S'^(t) and Jf , • • ■ ,Jn are mutually independent. 
Hence, Vw > and V6' = (6*1, • • • , 6^) > we have 



E [exp {-ujU{t) - 0V{t)} I Xn{t) =x]=E 



where = (-^("i)'--- ^J^) 



exp 



E 



exp 



Step 2 

Since S^(t) is the value, at time t, of a OLP with Levy measure H42|l we have 



E 



where 



exp {^.^(,}]^ exp {-^^04 
= /o'(l~exp{^y})A(dy) 

_ a fl l-cxp{-mi} 
- a" Jo 

= A(x)G„(l^) . 
On the other hand, using 143|l . Vy > 1 we have: 



Pi — >y 

X 



Hxy) ^ !_ 
A{x) y° 
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(45) 



Hence, combining the above together, we conclude that 

E,[eii-p{-LuU{t) -0V{t)} I X„{t)=x] = c^p{-A{x)Ga{uj)t} ■ L{0) , (44) 

where L{9) is the Laplace transform (at the point 9 ~ {9i, ■ • • , 6'„)) of the order 
statistics of {Yi, • • • , y„} which are IID Pareto(a)-distributed random variables: 
P(y >y) = l/y", y > 1. 

Step 3 

Using (T)), the probability density function of Xn+iit) is given by 

P(X„+i(t) e dx) 

= Pl{tX{dx))-Pr,{tJ^ X{ds)) 

= (tA(dx)).MMrexp{-tA(:.)} 

= i^exp{-tA(x)}A(x)"A(dx) . 
Combining 144|l and (|45|1 together we hence obtain 
E[exp{-wC/(t) - 9V{t)}] 

= E [exp {-iuUit) - eVit)} I X„(t) = a;] P (X„(i) e dx) 

= / • m , 

where 

j-n+l roo 

I = —- / exp {-A(x)GJLj)t} ■ exp {-tA(x)} A(xVX(dx) . 
nl Jo 

However, using the change of variable s = A{x) gives 

j-n+l roo 1 

/=-p/ exp{-t(l + G„(^))..}^"rf-s= , 
Jo (l + GQ(a;)) 

and hence we can conclude that 

E [exp {-LoUit) - 9Vm = (^^-i-^^ ^ . L{9) . (46) 

Equation gHl implies that {U{t): V{t)) = (^^ + h C„+i ; ^(i), • ' ' ,^(n)) 

where: 

(i) Ci7 ■ ■ ■ I Cn+iJ ^ij • • • J ^ are independent random variables; 

(ii) Ci; • • • iS.n+1 are C/a-distributed; 

(iii) Y(i), • • • , Y(„) are the order statistics of Yi, • • • ,Yn which are Pareto(a)- 
distributed: P(y > y) = l/y", y > 1. ■ 
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